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Abstract. The propagation of solitons in dipolar BEC in a trap potential with a 
barrier potential is investigated. The regimes of soliton transmission, reflection and 
splitting as a function of the ratio between the local and dipolar nonlocal interactions 
are analyzed analytically and numerically. Coherent splitting and fusion of the soliton 
by the defect is observed. The conditions for fusion of splitted solitons are found. In 
addition the dclocalization transition governed by the strength of the nonlocal dipolar 
interaction is presented. Predicted phenomena can be useful for the design of a matter 
wave splitter and interferometers using matter wave solitons. 
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1. Introduction 



The Bose-Einstein condensate (BEC) of chromium (^^Cr), where long-range dipolar 
interaction between atoms plays the dominant role, is a novel kind of nonlinear system 
becoming available to experiments [1]. Properties of dipole-dipole (DD) interactions, 
namely their long-range character and anisotropy, allow dipolar condensates to exhibit 
many unusual properties not found in BECs with just contact interactions [21 12]. In 
particular the existence of stable isotropic and anisotropic two-dimensional (2D) solitons 
has been predicted for such cold quantum gases [U [5] . Recently the bright solitons in 
quasi-lD dipolar BEC with competing local and nonlocal interactions have been studied 
in works [HllIllH]. 

The long-range dipolar interactions become dominant when the local part is detuned 
to zero by the Feshbach resonance (FR) techniques, as in the experiment on observation 
of Anderson localization in non-interacting cold quantum gases [9]. In this particular 
case the pure dipolar bright soliton can be observed. The propagation of such solitons 
under joint action of the trap and a barrier potential, including processes of crossing 
the barrier by the soliton during its oscillations in the trap are of the greatest interest 
to investigate. Such processes have a fundamental importance for the problems of 
entanglement of quantum solitons in cold dipolar gases. An interesting limit is the case 
of strong nonlocality when the dynamics of wavepackets become almost linear. 

This paper is devoted to the investigation of the properties of cold quantum gases 
in the presence of long-range dipolar interactions at the mean field level. First we study 
oscillations of the dipolar solitons in trap potential using the variational approach (VA) 
and the scattering theory. Second we consider the soliton transmission, reflection and 
splitting through a barrier placed in the center of the trap. Particular attention will be 
devoted to a process of coalescence of colliding wavepackets at the barrier which could 
be important for the design of beamsplitters and matter wave interferometers using 
matter wave solitons [10]. Finally, by means of numerical simulations of the original 
dynamical equation we verify the predictions of the VA and analyze the results beyond 
the analytical predictions. 



2. The model 



We consider the quasi- ID dipolar BEC loaded in a parabolic trap with a barrier at the 
center of the trap. The governing equation is the ID Gross-Pitaevskii equation (CPE) 
with a nonlocal interaction term [6l [11] : 

_ vI;(X,T) / rfe^(lX-el) |*(e,T)p = 0. (1) 

Here cuj, corresponds to the transverse trap frequency, /_l = ^JTl/ mu± , and d is the 
magnetic dipole moment oriented along the X-direction. The parameter e is connected 
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to the angle (p which the dipoles form with the longitudinal axis X. If dipoles rotate 
rapidly in the plane perpendicular to the axis X, the parameter e can vary from 1 
for dipoles oriented along the X-axis {(p = 0), to —1/2 in the case of perpendicular 
orientation to the X-axis ((/? = it/ 2) [6l [12]. The wave function is normalized to the 
number of atoms comprising the BEC, N' = |\E'(X)p(iX. The Hamiltonian for our 
model has the following form 



/oo 
(iX^*(X) 
-oo 



Vtr{X) + ha,u^\^iX)\' 



2m (9X2 



+ '-^Jd^^*{OR{X-0^iO 
Now we define dimensionless parameters: 



VI. (X). (2) 



ead 



t = Tu±, x = X/l±, g = - -, q=- -, = -y/2|a,o|^ 



1 1 ' 


g = 1 1 




\0'sO\ 



where = md'^/h'^ is the characteristic scale of the long-range dipolar interactions, and 
is the background value of the atomic scattering length. The dimensionless s-wave 
scattering length q is expressed in units of a^o due to the choice of the normalization of 
the wave function. Eq. ([T]) now can be written in the dimensionless form as follows: 

/+QO 
R{\x-^\) m,t)\'d^ = 0, (3) 
-oo 

where ilj{x,t) is the mean-field wave function of the condensate, q is the local contact 
interaction term, g{t) is the nonlinear coefficient responsible for long-range dipolar 
interactions, assumed to be time- dependent [12]. The external trapping potential Vtr{x) 
is 

Vtr{x) = VM + Vd{x), (4) 

VM = \u^x\ Vd{x) = Voe-^'/^'''\ (5) 

The frequency u is the longitudinal frequency of the trap, Vq and / are the amphtude 
and the width of the barrier. Potential Vd{x) for the case of a broad soliton. Is 3> I, can 
be approximated by a delta-barrier = a6{x) with a = V2ttIVq. The wave function 
ip{x) has normalization \ip{x)\'^dx = N. 

The following two forms for the kernels in the nonlocality term are possible 

R^(x) = (1 + 2x^) exp(x2) erfc(|x|) - 2Tr~^^^\x\, (6) 
R,{x)=xl{x' + xl)-'/\ (7) 

The former kernel corresponds to the dipolar BEC in a quasi-lD trap [6], while the 
latter one, which contains a cutoff parameter Xc also proposed for dipolar BEC in [8] 
(see also [13]) is more convenient for analytical treatment. 
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Using the matching conditions for Ri{x) 

/oo POO 
Ri{x)dx= / R2{x)dx, (8) 
-oo J —oo 

it can be found that Xc = tt"^/^. The comparison of profiles of the kernels ([6]) and 
at this choice of Xc shows very good fit. The cutoff parameter Xc has the meaning 
of an effective size of the dipole and the value is fixed by interpolating the function 
Ri{x) by R2{x) . Actually, it takes the value of the order of the transverse confinement 
length, which makes the model one-dimensional, and is the unit length in Eq. Q. In 
the limit x ^ Xc, where DD interaction effects dominate over the contact interaction 
effects, both response functions behave as ~ 1/x^. Thus, in the following we will use the 
kernel function R2{x) as it is more simpler for analytical treatments in the description 
of dipolar effects in BEC, using for example VA, where the integrals with R2{x) can be 
calculated in explicit form. 



3. Variational approach 

To describe the soliton propagation we employ the VA [T3]. The Lagrangian for Eq.([3]) 
has the following form: 

L = \{r^t - CO.) - - - VUxM' 

- r R{\x - mH^,t)\'dt (9) 

Z J —oo 

To derive equations for soliton parameters we use the Gaussian ansatz 

i;{x,t) = A{t) exp |-^^^^^ + - CW]' 

+ zk{t)[x-at)] + z(j){t)}, (10) 

where A, a, b, (, k, are the soliton amplitude, width, chirp, coordinate of the center of 
mass, wave vector and linear phase, respectively. Calculating the averaged Lagrangian 
L = J Ldx with this ansatz we obtain 

— = - + kCt-(ht r - b^a^ cjM — + C 

N 2 ^* ^* 4a2 2 2 \ 2 ^ y 

N 

[q + gF{a)]-V,G{a,U). (11) 



2V27ra 

where N = ^J^\A^^a =const and 



oo ^~OLt 



G(a, C, /) = , (13) 
/2 + aV/2 



In computing the integral the shift x— J-x — ^ ^^~C ^ind a transformation to 
the new variables z = l/2{x — C,), y = 1/2 (x + ,^) and dxd^ = 2dzdy have been used. 
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Considering the Euler-Lagrange equations 

d{L/N) d d{L/N) 



drji dt d{ri. 



•it 



we can derive the system of evolution equations for parameters rji = a,b,(,k,4' of the 
solution (fTOl) . Variations on (, k, b and a, respectively, give the following equations 

h = -uj\ - Vo^, (14) 

Ct = k, (15) 
at = 2ab, (16) 

- 1 25^ + «^ _ ^9 f M - i.^ - 1^^. (17) 



2a^ 2V27ia^ 2^f2^ada \a J 2 a da 

The equation for the phase is decoupled from the system, and we did not write it 
here. Finally, from ( fTTl) we obtain the system of equations for the soliton width and the 
center of mass 

where 

Ua = ^ + (g + gF) + + 2VoG, (19) 

Uc = ItoV + VoG. (20) 

Thus the evolution is described by the dynamics of two coupled nonlinear oscillators. 
From the first equation we can find the fixed point for the soliton width avA and from 
the second we can calculate the effective potential U^{avA) for the center of mass of 
the soliton. The description of the dipolar soliton dynamics by the system of two 
coupled nonlinear oscillators was applied successively, for example, in the work [T5] 
where soliton-soliton scattering in the dipolar BEC placed in unconnected layers was 
considered. This method was also used to observe oscillations of the solitons profile in the 
quasi-lD dipolar BEC in Ref.[7]. We also want to note that the Gaussian trial function 
gives good results for describing the stationary states and dynamics in quasi-lD and 
-2D geometries (see for example [T5l [7]). In some cases, for 2D geometries, when the 
purely dipolar BEC profile (for parameters unstable to collapse) has blood-cell forms, 
this ansatz failed, and it is necessary to choose the trial function as the sum of Gaussian 
profiles |T6l[T7]. In our case the failure of the ansatz can occur when the pulse is splitted 
into two or more parts by the defect. 
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4. Numerical results 

4.I. Stationary modes analysis 

We start with the case of competing nonlinearities {q ■ g < 0), namely we consider an 
attractive local nonlinearity, q < 1, and a repulsive DD interaction, g > 0, or vice versa. 
By using ansatz %p{x^t) = il)[x)e~^^^ and considering g{t) = g we get the stationary 
equation 

/+00 
^(k-el) \m\'d^ = 0. (21) 
-oo 

Let us first consider stationary solitons without external potential [uj = 0) and 
without defect (Vq = 0). In Fig{T] and [3] by solving numerically Eq. (l2T]) (by Newton 
iteration algorithm) we present the existence curves for families of the solutions with 
different combination of local and nonlocal terms. In FigJTJ^a) we start with the simple 
local case (g = — 1 and g = 0) for which the norm of the solution N goes to as the 
chemical potential fi goes to 0. By increasing gradually the coefficient g in the nonlocal 
term we pass through the critical value of g = gcr ~ 0.87 for which, at some fi = ficr, 
the existence curve starts to have a local minimum becoming bounded by critical value 
Ncr = N{ficr) which means that below N^r the solutions do not exist. For two branches 
in FiglT](a) corresponding to the purely local case {q = —1, g = 0), and with competing 
local and nonlocal terms, {q = —1, g = 0.8) we calculated numerically the width of the 
stationary solutions 

/oo 
x'\ij\'dx, (22) 
-oo 

and compared it with the results of the VA taken from the condition dUa/da = (see 
Eq. lfTO!) ). The results are shown in Fig{T]^b). For the purely local case {g = 0) one 
observes good agreement between the width calculated numerically and from the VA. 
However by increasing the repulsive nonlocal term the discrepancy between anum o-nd ay a 
starts to grow and for g > 0.8 results from the VA for small \fi\ does not agree with the 
width of the solution calculated numerically. 

In Fig|2] we have checked dynamically the stability of the solutions at the black 
and red points indicated in FigHt^a) where the derivative dN/dfi has different signs. 
According to Vakhitov-Kolokolov stability criterion these two solutions should have 
different stability. In FigJJt^a) the profiles of these two solutions are shown and in 
the panels (b), (c) their evolutions are presented. As it is expected the solution with 
dN/ dfi < is stable while the solution with dN/ dfi > is unstable, which is confirmed 
by direct numerical simulation in Figl2]J^b), (c). 

Similar analysis for the existence of the solutions with attractive nonlocal 
interaction {g < 0) and with the presence of the local repulsive interaction {q > 0) 
is presented in Figj3l Now, contrary to the previous case, the norm of the solutions is 
always unbounded which means that — ?■ as /i — ?■ 0. Also in this case the width of 
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Figure 1. In (a) existence curves N{^) for attractive local q = —1 and for repulsive 
nonlocal < .9 < 1 interactions. Solid and dashed lines correspond to the stable and 
unstable regions. In (b) the width calculated from the VA, ay a (solid lines), and from 
numerical solution of stationary problem, anum (points), are compared. 
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Figure 2. In (a) the profiles of the solution for q = —1, /.i = —0.1 and different values 
of coefficient of nonlocal term g = 0.8 (stable) and 5 = 1 (unstable) corresponding to 
the points in Fig[T] (a). In (b) and (c) the density plots of evolution of stable and 
unstable solutions, checked through direct dynamical simulation of Eq.([3]). 



the solutions calculated numerically from Eq.( l22|) and analytically from the VA are in 
very good agreement (see Figj3](b)). 




Figure 3. The same as in Fig (T] for attractive nonlocal interaction, g = —1, and for 
repulsive local interaction < q < 1. 
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5. Soliton dynamics in the trap and interaction with locahzed defect 

Now we consider the dynamics of the sohton with local and nonlocal interactions in 
the presence of the parabolic potential and the localized defect Va- Since we are 
interested in competing local and nonlocal nonlinearities (condition q ■ g < 0) in the 
following we will consider two possible combinations separately: i) g < 0, > 0; and ii) 
q>0,g<0. 

5.1. The case q < 0, g > 

Let us first concentrate on the case of attractive local (g = —1) and repulsive nonlocal 
{g = 1) interactions when the existence curve is bounded by the critical norm Ncr, which 
occurs above g ~ 0.87. In FigJUthe initial profile of the unperturbed stable soliton taken 
from the existence curve in Fig|T]at /i = — 1 and shifted to the position xq = 4, as well 
as the corresponding external parabolic potential with the defect are shown. The shift 
from the center is used to get oscillations of the soliton in the defect-free parabolic 
potential. 



- l¥l' 
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Figure 4. Initial profile of the soliton \ip{xQ — A,t — 0)| (solid line) in the parabolic 
trap with the defect Vd (dashed line). The parameters are: ii ^ —1, q = —1, g = I, 
Vo = 1, e^O.l, Lo = 0.2. 

As soon as the soliton is shifted from the minimum of the parabolic potential it 
starts to accelerate to the center and the soliton velocity at the center, Vc, will depend 
on the strength of the parabolic potential u and magnitude of the shift Xq according to 
Vc = ojxq. In this way by changing the strength of the trap with fixed Xq we can control 
the soliton velocity in the process of the interaction of the soliton with the defect. In 
the defect-free parabolic trap the frequency of oscillations of the soliton coincides with 
the trap frequency {((t) = cos(a;t)). By fixing the amplitude Vq and the width £ of 
the defect and letting the soliton collide with the defect we found three characteristic 
regimes of the soliton-defect interaction depending on the incoming velocity. In FigJS] 
we present two regimes. The first one corresponds to a relatively small soliton velocity 
(weak trap) when the soliton is refiected from the defect and becomes "closed" in the 
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semi-space (x > 0) of the parabolic potential (see FigJSJ^a)). By increasing the soliton 
velocity by changing the strength of the parabolic potential we found another limiting 
case when the soliton has enough kinetic energy to pass through the defect (see FigjSj^b)). 
To visualize the corresponding evolutions we calculated the norms in the right part and 
in the left part from the defect as Nr = \'ijj\'^dx and Ni = N — Nr. The results are 
shown in Figl5t^c),(d) which confirm the above mentioned soliton behavior. 

We also compared direct numerical calculation of the soliton dynamics with the 
results obtained from the VA. In Fig|5]^a),(b) by dashes lines we present the trajectories 
of the oscillations of the center of mass of the soliton calculated from Eqs. (fT9|) . ( 120|) . 



ro=0.05 



(a) 





w=0.5 ' (b) 








-n 



2 : 
ok 



100 



Figure 5. In (a), (b) the density plot of the osciUations of the sohton in parabohc 
potential in the presence of the defect. The strength of the parabolic potential in (a) 
w = 0.05 and in (b) to = 0.5. Initial profile of the soliton is shifted to xq = 4. The 
dashed black lines correspond to the trajectories of the soliton calculated from the VA 
(fT9| . (|20l). Other parameters are: fi = -l,q = -1, g^l,Vo = l,e = 0.1. The dotted 
line in (a), (b) shows the position of the defect. In (c) and (d) the evolutions of the 
norms to the right, Nr (sohd black), and to the left, A^; (dashed red), from the defect 
are shown corresponding to the cases (a) and (b). respectively. 



A more interesting effect of the coherent splitting of the soliton by the defect is 
observed for intermediate values of u. This is the case when the VA fails and only 
numerical calculation will be presented. We found (see Figj6]) that taking the strength 
of the parabolic potential around u ^ 0.2 one can observe splitting of the soliton into 
two quasi-symmetrical parts. This can be seen from the evolution of the norms Nr and 
Ni in Figj6]^b), (c). It should be stressed here that after splitting the soliton lost its 
solitonic identity as soon as in this case the existence curve has bounded norm Ncr and 
the norm of the each parts is below this critical value (A'^; < N^r and A^,. < Ncr). To 
confirm this behaviour we switch off the parabolic trap and take a soliton with initial 
velocity Vc ~ 0.8 far from the defect. As it is shown in FigjTl after interaction with the 
defect the soliton splits into two packets which transform continuously into linear waves. 
In the presence of the trap these linear packets do not escape but they are reflected by 
the parabolic potential and return to the center where again produce the initial soliton 
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Figure 6. The same as in Figl5]witli u = 0.2 in (a) and ui = 0.203 in (b). 



("fusion") in the left or right part of the defect and continue this process through several 
periods. In Figj6]^a) one observes that after splitting and returning to the center the 
soliton continues to move towards the negative x passing completely through the barrier 
while in[6t^b) soliton after fusion is reflected from barrier. It should be stressed that very 
tiny changes in the initial conditions could affect the condition for the soliton splitting 
and the dynamics of the soliton after fusion (transmission or reflection scenarios). 

In this case the VA can be used to find the condition for soliton splitting. As follows 
from Eqs.flTTl). (|T9l) . when the kinetic energy of the effective particle Ekm = x\uj^ jl is 
less then the barrier effective potential height E^ar = VoG(avA5 ^5 C = 0), a full reflection 
of the soliton occurs. In the opposite case we have full transmission (compare FigsJS]). 
Therefore we can assume that when E'^j^ = E^y^r^ the partial reflection/transmission 
should take place. Considering this condition as the soliton splitting condition, we 
obtain: 

xo = ^v'2K,G(ayA,^,C = 0). (23) 

By comparing the numerical results with the analytical ones taken from (123 p one observes 
excellent agreement (see Figj8]). 

For smaller strength of the repulsive nonlocal nonlinearity, when the norm of the 
solution is unbounded, one observes a rather different picture of splitting of the soliton 
at the defect (see FigsjQ]) 

Comparing Figj6] and Figl9] one can conclude that in the case g = 1 the soliton 
transforms into two quasi-identical wave packets in Figj6]|^a),(b) (one also observes this 
effect looking at A*";^^ when along a half period A"; ^ Nr) while in the case g = there is 
no splitting of the soliton. Instead of that it transforms into the defect mode localized 
at the defect. 

This can be verified in FigJTO]by comparing the existence curves for the case without 
defect with curves calculated in the presence of the defect. As one can see in FigJTOlby 
switching the defect on the existence curves go upper and what is essential is that the 
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Figure 7. Scattering of the soliton with initial velocity Vc = 0.8 placed at a;o = 4 in 
the absence of the parabolic trap (oj = 0). Other parameters are the same as in FiglS] 




Figure 8. Comparison of numerical (line with points) and analytical (dashed line) 
results using Eg. ([^5]) for the critical coordinate at which one can observe splitting 
effect {Ni « Nr). Parameters are: w = 0.2, ^ = -1, g = -1, .g = 1, £ = 0.1. The 
width of the soliton calculated from the VA is ay^ w 0.38. 
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Figure 9. Dynamics of the soliton with q = —1 and g ~ in the parabolic trap with 
the strengths uj = 0.145 (a),(c) and uj = 0.15 (b),(d). 
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critical value of g at which the existence curves become bounded also decreases with 
presence of the defect. As an example the existence curve for g = 0.8 in the defect-free 
case, Vo = 0, is unbounded while in the presence of the defect with Vq = 1, the existence 
curve becomes bounded. Considering the case g = the soliton at yU. = — 1 and Vq = 
has the same number of particles as the defect mode at = —0.68 and Vq = 1 (see the 
lower horizontal dotted line in the left panel of Fig JTOj) . To check this we calculated the 
defect mode for fi = —0.68 in the presence of the defect and compared it with the defect 
mode obtained by direct dynamical calculation where the initial soliton for /i = — 1 and 
Vq = has the same number of particles as the defect mode. In the right panel of FigJTOl 
we compare these two defect modes (time t = 8.37 corresponds to the instant when 
A^^. ^ A*"; during interaction of the soliton with the defect). 

In the case g = 0.8 the situation is different. The number of particles in the soliton 
at /i = — 1 and Vq = is below the critical number of particles, Ncr, needed to generate 
the defect mode in the presence of the defect (there is no intersection of the horizontal 
dotted with the upper dashed line). In this case one can observe splitting of the soliton 
into two linear wave packets after interaction with the defect. 
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Figure 10. Left panel: Existenee eurves without defeet Vq = (solid lines) and 
with defect Vq ~ 1 (dashed lines) for the cases 5 = and g = 0.8. Right panel: 
The profiles of the defect modes with the same number of particles calculated form 
stationary problem (black solid) at the black point in the left panel with fi = —0.68 
and from dynamical equation (dashed red) where initial solution was taken at fj, — —1 
for Vo = 0. 



Localization- delocalization transition governed by the nonlocal term. Let us consider a 
linear variation in time of the strength of the nonlinear nonlocal coefficient of the form 

g{t)=gf + ig^-gf)\l-2t/tf\ (24) 

where gi = g(t = 0) and gf = g{t = tf/2) correspond to the initial strength of the 
nonlocal coefficient and its value at the turning point, t = tf/2, with the return to 
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the initial value at t = tf. Tuning the dipolar interactions in quantum gases can 
be achieved using for example time dependent control of the anisotropy of dipolar 
interactions suggested in [12]. As it was shown in Section 4.1, in the presence of 
the attractive local and the repulsive nonlocal interaction (see Fig{T]) one can have a 
transition between unbounded and bounded cases of the existence curves. Using a time 
dependent nonlocal coefficient g{t) in Fig. [11] we present evolutions of the profile of 
the soliton density with \gf\ < \gcr\ and \gf\ > \gcr\ (here gcr corresponds to the case 
when the norm of the initial soliton coincides with the norm at the global minimum of 
the existence curve N^r)- As one can see in the former case there is no delocalization 
transition and the solution remains localized, while in the latter case at some instant 
the norm of the solution becomes smaller then the critical norm N^r and the norm of 
the solution any more pretends to the existence curve and therefore the solution decays 
into linear waves. Returning to the initial value g = gi in the first case the solution 
reconstructs its initial form while in the second case it remains delocalized. 



8 
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Figure 11. Localization-delocalization transition managed by nonlocality. In (a) the 
amplitude of the density in the center \ip{x = 0)p for two different values of gj: black 
solid line 5/ = 0.85 and red dashed line 5/ = 0.95. In (b) and (d) the density plot of 
the corresponding dynamics. Other parameters: /x — —1, q = —1, gi = 0.8, w = 0, 
^0 = 0. 

5.2. The case g < 0, q > 

Now let us consider the opposite situation when the solitonic structure is supported by 
an attractive nonlocal interaction g < and the local interaction is repulsive, g > 0. 

In FigUl] the corresponding initial profile of the nonlocal soliton in the parabolic 
trap with defect is shown. As in the previous case the shift from the center of the trap is 
needed to observe oscillations of the solitons and eventual interaction with defect placed 
in the center of the trap. 
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Figure 12. Initial profile of the soliton |^/'(xo = 4,t = 0)| in the parabolic trap K; 
with the defect Vd- The parameters are: /i = —1, q = 1, g = —1, Vq = I, £ ^ 0.1, 
uj = 0.15. 



Two simple scenarios of the interaction of the sohton with the defect in the case of 
a strong parabohc trap are shown in Fig JTST a). (c). In the first case the incoming sohton 
has enough kinetic energy to go through the defect. In this interaction the soliton lost 
part of the energy and therefore became locked by the defect in the semi-space x < 0. 
For a stronger trap the soliton has enough kinetic energy to go through the barrier 
several times (see Fig JTST b). (d)). 




Figure 13. Dynamics of the dipolar soliton with q = 0.5 and g = —1 in the parabolic 
trap with presence of the defect Vb = 1, ^ = 0.1. In (a), (c) the evolution of the density 
plots are shown for uj = 0.14 in (a) and a; = 15 in (b) with corresponding dynamics of 
the norms Nr and Ni. 



A similar picture of the scattering of the pure dipolar soliton can be observed by 
decreasing the strength of the local interaction q to zero (see Fig JT^ . As one can see 
the pure dipolar soliton does not present very robust dynamics, what can be explained 
by the weak force of the nonlocal interaction. 
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Figure 14. The same as in FiglT3]for the case of the pure dipolar sohton with q ~ 
and g = —1. The strength of the parabolic trap in (a), (c) is a; = 0.14 and in (b),(d) 
is u! — 0.15. 

6. Conclusion 

In conclusion, we have investigated the dynamics of bright solitons in a dipolar 
condensate loaded into a parabolic trap with a barrier potential at the center. Using a 
variational approach and scattering theory, we have studied the reflection, transmission 
and splitting of bright solitons on the defect in the presence of the trap potential. We 
explore the different sets of parameters changing the relative strength of the local and 
nonlocal (dipolar) interactions. The case when the local nonlinearity dominates is close 
to the settings investigated early [HI [191 120] and well understood now. The intermediate 
cases when q ^ g and the pure dipolar nonlinearity g = are considered here. 

We show that coherent splitting of the dipolar bright soliton exists. The splitting 
parts are reflected by the trap and recombine into a single soliton, performing a few 
oscillations under the trap potential. The condition for soliton splitting is derived, 
which is in excellent agreement with numerical simulations of the full nonlocal Gross- 
Pitaevskii equation. 

Localization-delocalization phenomena governed by the variation of dipolar 
interactions in time has been studied for the cases \gf \ < \gcr\ and \gf\ > \gcr\- The 
reconstruction of the soliton to its initial form is observed for the former case, while in 
the latter case the soliton transforms into the linear waves. 

The obtained results can be interesting for the study of soliton dynamics in the 
case of very small atomic scattering length when the dipolar interactions effects start to 
play a dominant role. Several experiments were performed recently in ^Li with a very 
small atomic scattering length tuned via the broad 1 1 , 1 > Feshbach resonance with the 
number of atoms in the soliton of the order of ~ 2 ■ 10^ [101 [21]. The barrier potential 
was generated by a near-resonant cylindrically focused laser beam. It can be useful 
for the design of matter-wave beamsplitters and matter-wave interferometers using the 
fusion of the solitons. 
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The related problem of quantum scattering of solitons by the potential well for 
BEC with local interactions has been studied in [221 [231 |2l]. For the case of the dipolar 
solitons this problem will be considered separately. 
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